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Abstract 



The inclusive energy spectra of light charged particles, such as, a, p, d and t, 
evaporated from hot ^^P nucleus at excitation energy E* ~ 60 MeV, have been 
measured at various angles. The compound nucleus ^^P has been populated 
using two different entrance channel configurations; i.e., ''Li (47 MeV) + ^^Mg 
and ^^F (96 MeV) + ^^C reactions, leading to same excitation energy of the 
compound system. It has been observed that the spectra obtained in the ''Li 
(47 MeV) + ^^Mg reaction follow the standard statistical model prediction with 
a spherical configuration of the compound nucleus. But, the spectra obtained 
in the ^^F (96 MeV) + ^^C reaction deviate from similar predictions of the sta- 
tistical model both on higher as well as on lower energy sides. Considerable 
deformation was required to be incorporated in the calculation in order to re- 
produce the measured energy spectra in this case. Dynamical trajectory model 
calculations were not found to play any significant role in explaining the differ- 
ences in behaviour between the two cases under study. The observed discrepancy 
has been attributed to the difference in the angular momentum distributions of 
the compound nuclei formed in the two reactions. 
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I. INTRODUCTION 



In recent years, there has been a lot of interest to study the statistical properties of hot 
rotating nuclei produced through heavy ion fusion reactions. High excitation energy favours 
the nucleus to deexcite through statistical emission of several light particles followed by 7 
decay. The deexcitation process is governed by the available phase spaces, i.e., the nuclear level 
densities of both the parent and daughter nuclei, as a function of excitation energy, angular 
momentum and deformation. The evaporated light charged particle spectra in conjunction with 
the statistical model calculations can be used to extract the statistical properties of the hot, 
rotating nuclei. 

However, recent studies reveal that, at high excitation energy and angular momentum, the 
energy spectra of evaporated light charged particles, are no longer consistent with the standard 
statistical model predictions. A large number of experiments have been performed in last few 
years to study this anomaly over a wide range of compound nuclear masses {^cn) in the range 
~60-170, for example [0-0, and in most cases (with the exception of [Q), it has been found 
that the statistical model fails to explain the observed light charged particle energy spectra. 
Choudhury et al have measured the a-particle energy spectra in the reaction 214 MeV ^^S on 
^''^Al and observed that the measured spectra deviate considerably from the standard statistical 
model predictions. La Rana et al [§] have measured a-particle spectra in the reaction 190 MeV 
^°Ar on ^''Al and have seen similar discrepancy. The authors argued that in the statistical 
model calculation, the emission barrier is considered as the fusion barrier of the two partners in 
the exit channel assuming their spherical configuration. It may be a good approximation at low 
excitation energy and angular momentum. But, when the nucleus is at high excitation energy 
and rotating with high spin, it is likely to be deformed in shape. Hence, the emission barrier will 
be lower if one considers that the a-particles are emitted preferentially along the larger axis. On 
the contrary, it has been claimed that, only deformation, unless it is unrealistically large, 
cannot affect the a-particle emission barrier much; however, the anomaly may be well explained 
by incorporating spin dependent level density in the standard statistical model prescription to 
modify the yrast line keeping the emission barrier unchanged. Recently, Charity et al ^ have 
measured light charged particle spectra emitted in the reactions ^^Ni on ^''^Mo and ^^O on ^^^Sm. 
They have shown that the measured spectra deviate from the statistical model predictions at 
incident energies above ~5 MeV/nucleon, and that different set of parameters are required to 
reproduce the data for different types of light charged particles, such as, a, proton, deuteron and 
triton. Moreover, magnitude of the deviation has been found to depend on the excitation energy 
as well as the entrance channel. Intuitively, the entrance channel dependence may affect the 
shapes of the light charged particle spectra if the equilibrating system remains deformed over a 
time scale comparable to the mean lifetime of charged particle emission. Recent measurements 
of a-particle energy spectra for the reactions 83, 97 MeV ^^C on ^^Sc |Tl| and llOMeV ^^O 
on ^'^Fe |]10| have also indicated that the role of entrance channel dynamics should be carefully 
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looked into. 

It is thus apparent that there exists some discrepancy between the experimentally measured 
spectra and the corresponding predictions of standard statistical model, which is prominent for 
higher spin and for heavier systems (A>60). However, due to lack of experimental data situation 
is not clear for the lighter systems at relatively lower spin. This prompted us to explore the 
scenario for light system, such as, ^^P, in the incident energy range below 5 MeV/A. In the 
present work, we have populated the compound nucleus '^^P at excitation energy ~ 60 MeV 
through different entrance channels, i.e., 96 MeV ^^F on ^^C and 47 MeV ^Li on ^^Mg and 
have measured the energy and angular distributions of all light charged particles, (i.e., a, p, 
d and t). Since the two reaction channels populate the same compound system with different 
angular momentum distributions (Icr ^ 21 h and 16 h, respectively), it offers an opportunity to 
study the angular momentum dependence in the evaporative decay cascade. Apart from that, 
entrance channel configurations may also affect the evaporative decay if the nucleus lives in 
the dinuclear form for a long time and evaporates during this period. But our previous studies 
on the intermediate mass fragment emission did not reflect the existence of any such long- 
lived dinuclear configuration |[13|Jl^ of the ^^P compound nucleus. Earlier, we have measured 
also a-particle spectra in coincidence with the evaporation residues in the 96 MeV ^^F on ^^C 
reaction and have observed that the measured spectra deviate from the statistical model 
predictions both at low as well as at high energy sides. Assuming a deformed configuration of 
the compound nucleus, the full energy spectra could be explained satisfactorily. Interestingly, 
it has been observed in earlier works that, different amount of deformations are required to 
explain different sets of light charged particle data. It directly points to the inadequacy of the 
standard statistical model to explain the behaviour of light charged particle emission spectra 
with increasing spin and excitation energy. In the present work, we would like to explore 
systematically the scenario of light charged particle emission further for a lighter system at 
moderate spin and excitation energy. 

The paper has been arranged as follows. The experimental details alongwith the experi- 
mental results are described in Sec. |T[ The predictions of the statistical model for ^^P nucleus 



have been discussed in Sec. |T|. Finally the summary and conclusions are presented in Sec. IV 



II. EXPERIMENTS AND RESULTS 

The experiments were performed at the Bhabha Atomic Research Centre - Tata Institute 
of Fundamental Research 14 UD pelletron accelerator laboratory, Mumbai, with 96 MeV ^^F 
beam on 125//g/cm^ self-supporting ^^C target and 47 MeV ^Li beam on 275 yug/cm^ self- 
supporting ""^Mg target. The beam current was typically 20-80 nA. The light charged particles 
have been detected in three solid state telescopes. The telescopes (Ti, T2, T3) were of two 
elements, consisting of lOO/xm, 45/im, 40/im AE Si(SB) and 5 mm, 5 mm, 2 mm Si(Li) E 
detectors respectively. Typical solid angles were 1.5 msr, same for all the three telescopes. 
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Analog signals from the detectors were processed using the standard electronics before being 
fed to the computer for on-line data acquisition. 

The telescopes were calibrated with the radioactive ^^^Th a-source. Typical energy reso- 
lutions (FWHM) obtained were 1.12% for Ti and 1.13% for T2 and T3 respectively. The low 
energy cutoffs thus obtained were typically 3 MeV for proton, 5 MeV for deuteron, 5 MeV 
for triton and 11 MeV for a in Ti. In T2 and T3 cutoffs were 1 MeV for proton, 3 MeV for 
deuteron, 3 MeV for triton and 7 MeV for a-particles. 

Inclusive energy distributions for a-particles, proton, deuteron and triton have been mea- 
sured in an angular range of 25° - 140° in case of ^Li (47 MeV) -|- ""*Mg reaction and in a 
range of 5°-60° in case of ^^F (96 MeV) + ^^C reaction. In the centre of mass (c. m.) frame it 
covers up to ~ 140° in both the cases. Figs. show the measured energy spectra for the light 
charged particles at different laboratory angles respectively for the two reactions. The shape 
of the energy spectra is similar for all the measured angles indicating the fact that the thermal 
equilibrium has been achieved in both the cases before the particle emission. 

The average velocities of the observed light charged particles have been plotted in the parallel 
(f ll) and perpendicular {v±) velocity plane ( Figs. for the "^Li (47 MeV) + ^''Mg and ^^F (96 
MeV) + ^^C reactions. It has been observed, that for all particles, the average velocities fall on 
a circle around the compound nucleus velocities (fen)- It implies that the average velocities (as 
well as the average kinetic energies) of the light charged particles are independent of the centre 
of mass emission angles in both the reactions. It means that the energy relaxation is complete 
and the light charged particles are emitted from a fully equilibrated source moving with the 
velocity Vcn- The deviations in average velocity of the light charged particles at backward angles 
i^iab > 40°) are due to large energy cutoffs in the energy spectra. 



III. DISCUSSION 



The experimental data have been compared with the predictions of standard statistical 
model calculation using the code CASCADE |l^. Figs. H - H show the centre of mass 



energy spectra for the various light charged particles, such as a, p, d and t, observed in the 
^Li (47 MeV) + ^"^Mg reaction. The solid lines represent the predictions of the statistical 
model code CASCADE. The optical model parameters used for calculating the transmission 
coefficients were taken from [0 for p, d, t, and from [0 for a-particles. The critical angular 



momentum for fusion, Icr, is calculated using Bass model as 16h |T^. Sharp cut-off angular 
momentum distribution has been assumed. The radius parameter was taken to be 1.29 fm 
following Piihlhofer et al |T^. The angular momentum dependent deformation parameters {61 
and 62 as described later) were set equal to zero. From the figures, it is observed that the 
theoretical predictions match well with the experimental results, indicating that all the light 
charged particles are evaporated from a nearly spherical ^^P compound nucleus . 
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However, in case of ^^F (96 MeV) + ^^C reaction, the statistical model calculation using the 
code CASCADE, with the standard parameter set, fails (the dashed lines) to predict the shape 
of the experimental spectra. Figs. || - |l^ show the observed (filled circles) and calculated (solid 



and dashed lines) centre of mass energy spectra for the various light charged particles emitted 
in the ^^F (96 MeV) + ^^C reaction. The figures indicate deviation of the experimental data 
from the standard statistical model predictions (dashed lines) at low as well as at high energy 
sides of the spectra. 

According to the CASCADE calculations, the lower energy part of the spectra is controlled 
by the transmission coefficients. For the calculation of transmission coefficients, we have used 
the same sets of optical model parameters as in the case of ^Li (47 MeV) + ^^Mg reaction. 



{i.e., the optical model parameters from for alpha particles, and from |]T6[ for proton, 
deuteron and triton). A deformed configuration of the compound nucleus has been assumed in 
the calculation. This is realized by varying only the radius parameter ro. The optimum value of 
ro is found to be ~ 1.56 fm which is consistent with our earlier work 0. The increased radius 
parameter lowers the emission barrier and hence enhances the low energy yield by increasing 
the transmission coefficients. 

The higher energy part of the spectra is controlled by the level density p{E,j), which for a 
given angular momentum j and excitation energy E is defined as |]15[ , 



p^E,j) = '^l±la'/\J^f^- 1 -exp[2[a{E - A - EA^ (1) 

where a is the level density parameter, which is taken to be A/8 for the present calculation, 
T is the thermodynamic temperature, A is the pairing correction, and Ej is the rotational 
energy which can be written in terms of the effective moment of inertia jT^jj as 

i?, = -^j(j + l)- (2a) 

In hot rotating nuclei, formed in heavy ion reactions, the effective moment of inertia may 
be spin dependent and may be written as, 

J'e// = JoX(l+5i/ + 52/), (2b) 

where, the rigid body moment of inertia, j7o, is given by, 

Jo = \a''I\1. (2c) 
5 

Here, A is the mass number and ro is the radius parameter. Non zero values of the parameters 
5i, ^2 introduce spin dependence in the effective moment of inertia (Eq. [2bD . The relation (Eq. 
2a) defines a region in the energy, angular momentum plane (E-j plane) of allowed levels and 
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a region where no levels occur. The region is bounded by the yrast line defined by the lowest 
energy state that occurs for a given angular momentum. It has been observed in the previous 
experiments |^J^ that the yrast line has to be modified to explain the experimental data. 
Different attempts have been pursued using different tools to modify the yrast line, such as, 
radius parameter , deformation parameters Si , 62 and the angular momentum j . 

However, we have observed that the variation of a single input parameter, Tq, in the code 
CASCADE, is sufficient to take care of the low energy as well as the high energy side of the 
spectra 0. So, in the present calculation, we varied the parameter To only, and kept the spin 
dependent parameters 61, 62 in Eq. (pBD equal to zero, to reproduce the experimental spectra. 
It is clear from Figs. ^ - that the calculated spectra (solid lines) agree quite well with the 
experimental spectra both in the low energy as well as in the high energy regions. 

For the sake of completeness, we also attempted to interpret the data using the other 
alternative; i.e., by introducing angular momentum dependence in the effective moment of 
inertia through the angular momentum dependent deformation parameters, i.e., 61, 62 of Eq. |2b[ 
The best fit values 61 and S2 obtained from fitting the data were = 2.8 x 10~^ and S2 = 
2.5 X 10^^. A representative theoretical prediction with the above best-fit parameters have 



been shown in Fig. |T3| along with the respective experimental data for the a-particles emitted 
in ^^F (96 MeV) + ^^C reaction. It is clearly seen that with the above values of 61 and 62, 
the predicted spectra match well with the experimental ones. It is further evident (comapring 
Fig. 13 with Fig. |^) that both the prescriptions (either changing Tq only or varying 61 and ^2) 
yield nearly equivalent result, at least for the system considered here. 

In order to probe more deeply into the matter, we have estimated an average radius param- 
eter (for non zero values of 61 and ^2) by averaging over the angular momentum distribution 
(sharp cut-off triangular distribution has been assumed, upto the critical angular momentum 
21h) in the following way. 



,2 _. n-a+Sij'+s2j%2j + i) 



Here, ro is taken as 1.29 fm. jcr is the critical angular momentum of the system. Si and 
82 values are as stated above. The effective radius parameter, re//, thus obtained is 1.65 fm 
which is within 10% of the radius parameter we have used in the statistical model calculation 
to reproduce the experimental data. Therefore, it can be said that the variation of single 
parameter ro to fit the data is nearly equivalent to the spin dependent level density prescription 
as described above, at least for light compound systems involving lower spins. 

For such light compound systems like ^^P at low spins, the present method has a distinct 
advantage. Here, both lower and higher energy side of the spectra could be fitted by changing 
a single parameter r^ once for all. The increase of ro causes lowering of the barrier which 
enhances the transmission coefficient, thereby correcting the lower energy side of the spectra; 
at the same time, it also causes a reduction of rotational energy thereby changing the level 
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density in right direction to explain the high energy part of the spectra. On the contrary, in 
the standard spin dependent level density prescription, one has to take care of the low energy 
side of the spectra separately by introducing a reduced barrier. However it may be noted 
that, as the spin dependent deformation becomes more significant at higher spins, the present 
method may not lead to realistic values of ro for higher spins; in such cases one may have to 
consider the spin dependent deformation only. On the other hand, for lower spin values, (for 
example, the compound nucleus ^^P formed in the reaction ^Li (47 MeV) + ^^Mg with Icr ~ 16 
h) the change in ro may not be quite appreciable. This intuitively explains why the standard 
statistical model predictions have been in good agreement with the data in this case. 

We have also tried to estimate the effect of entrance channel dynamics on light charged 
particle emission using the code HICOL . Fig. |T^ shows the calculated temparature of the 
fused nuclei as a function of time for various angular momentum values covering the range of 
interest for the two reactions under study. From Fig. |T^, it is observed that for both the cases, 
thermal equlibration of the system is achieved, on the average, within 8xl0~^^ second after 
the zero time. The zero time is defined as the time when the participating nuclei begin to feel 
the nuclear forces and deviate from their earlier coulomb trajectories. It is also observed that 
the temperature decreases with increasing angular momentum. Decay time has been estimated 
to be ~12xl0~^^ second, using the code PACE2, which is comparable to the formation time 
of the nucleus as obtained from HICOL. This implies that the effect of entrance channel on 
particle emission should be similar in both the reactions. However, particle spectra observed 
in the two reactions are different, implying that the above difference may not be due to the 
entrance channel dynamics. 



IV. SUMMARY AND CONCLUSION 



We have measured the inclusive light charged particle energy spectra in an angular range 
of 5°-60° and 25°-140° respectively in the ^^F (96 MeV) + and ^Li (47 MeV) + "''*Mg 
reactions. It has been observed from the average velocity plots in the v\\ vs. v± plane that the 
average velocities of all light charged particles in both the reactions stated above fall on the cir- 
cles centered around their respective compound nuclear velocities, implying that these particles 
are emitted from equilibrated compound nuclear sources. The data in the two reactions have 
been compared with the statistical model calculations using the code CASCADE. It has been 
observed that the light charged particle spectra measured in the ^Li (47 MeV) + ^^Mg reaction 
can be well explained in the framework of standard statistical model calculation, whereas the 
light charged particle spectra obtained in the ^^F (96 MeV) + ^^C reaction deviate significantly 
from the similar calculations. The formation times as estimated by the semiclassical calcula- 
tion using the code HICOL and the decay times as estimated by PACE2 are comparable in 
both the cases. Therefore, the discrepancy in the light charged particle spectra obtained in 
the two reactions may not be due to the entrance channel effect as the estimated formation 
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and decay times of the compound nucleus in both the cases are of the same order. Rather, it 
may be attributed to the difference in the angular momentum distributions of the compound 
nuclei formed in the two reactions considered here. A satisfactory description of the measured 
energy spectra in the case of ^^F (96 MeV) + ^^C has been achieved by invoking a deformed 
configuration of the compound nucleus through the modification of the radius parameter Tq 
from 1.29 fm to 1.56 fm in the statistical model calculation. The data could also be reproduced 
equally well by incorporating spin dependent level density formalism and the average radius 
parameter obtained using this prescription is 1.65 fm, which is within 10% of the radius pa- 
rameter we have used to reproduce the experimental data. However, it may be mentioned here 
that the present simphfied approach may not lead to realistic values of Tq for larger values of 
compound nucleus spin, which may limit the use of this simple approach for heavier compound 
systems at higher spins. In conclusion, the present work demonstrates that a modified radius 
parameter (ro) alone may be sufficient to reproduce the deformations for lighter systems (A ~ 
30). However, it will be interesting to explore the validity of this simple approach in explaining 
the light charged particle data for heavier systems at higher excitations. 
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FIGURES 




FIG. 1. inclusive (filled circles) light charged particle energy spectra obtained in the ^Li (47 MeV) 
+ ^^Mg reaction at different laboratory angles. For each curve cross-sections at highest observed angle 
is multiplied by 10° and the multiplication factor subsequently increases by a factor of 10 for each of 
the successive curve as the angle increases. 
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FIG. 2. Same as fig. 1 for tfie ^^F (96 MeV) + reaction. 
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FIG. 3. Average velocities of the light charged particles, obtained in the ^Li (47 MeV) + ^^Mg 
reaction, as a function of velocities parallel {v\f^ and perpendicular (v^ to the beam direction. The 
arrow indicates the velocity of the compound nucleus. 
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FIG. 5. inclusive (filled circles) a-particle energy spectra, obtained in the ^Li (47 MeV) + ^^Mg 
reaction, in centre-of-mass frame measured at different laboratory angles. Solid lines (normalized at 
the peak of the measured data) are the prediction of CASCADE calculations with the standard value 
of the radius parameter, tq. (see text). 
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FIG. 6. Same as fig. 5 for protons. 
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FIG. 7. Same as fig. 5 for deuterons. 
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FIG. 8. Same as fig. 5 for tritons. 
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FIG. 9. inclusive (filled circles) a-particle energy spectra, obtained in the ^^F (96 MeV) + ^^C 
reaction, in centre-of-mass frame measured at diflFerent laboratory angles. Curves are the predictions 
of CASCADE calculations with the standard (dashed curves) and the modified (solid curves) values 
of the radius parameter, ro (see text). Curves are normalized at the peak of the measured data. 



18 




FIG. 10. Same as fig. 9 for protons. 
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FIG. 11. Same as fig. 9 for deuterons. 
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FIG. 12. Same as fig. 9 for tritons. 
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FIG. 13. Inclusive cm. a-particle energy spectra at different laboratory angles, obtained in the 
^^F (96 MeV) + ^^C reaction. Solid lines (normalized at the peak of the measured data) are the 
predictions of CASCADE with Si = 2.8 x 10"^ and S2 = 2.5 x 10"'^ (see text). 
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FIG. 14. Time evolution of the temperature (T) parameter as obtained from HICOL calculations 
for the reactions studied. 
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